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Resumo 


Escoamentos viscosos nao sao, em geral, discutidos de maneira detalhada em 
cursos de fi'sica geral e basica. Isso se deve, em parte, ao fato da equagao 
de Navier-Stokes admitir solugao analitica apenas para alguns poucos cases 
restritos, enquanto problemas mais sofisticados so podem ser resolvidos por 
meio de metodos numericos. Neste trabalho, apresentamos um tunel de vento 
simulado, i.e., apresentamos um conjunto de programas capazes de resolver a 
equagao de Navier-Stokes para um objeto de formate arbitrario inserido em 
um tiinel de vento. O tunel nos possibilita visualizar a formagao de vortices 
atras do objeto, os contiecidos vortices de von Karman. E possivel tambem 
calcular a forga de arrasto sobre o objeto. Acreditamos que este tiinel de vento 
numerico possa subsidiar o professor e permitir que uma discussao mais elabo- 
rada de escoamentos viscosos seja realizada. As potencialidades do tiinel sao 
exemplificadas por meio do estudo do comportamento da forga de arrasto so¬ 
bre um modelo simplificado de asa cujo angulo de ataque pode ser controlado. 
Um link para o download dos programas que compoem o tiinel e apresentado 
ao final do texto. 


* paulo. victor @ifrj. edii. br 
Vmco@if.uff.br 
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Abstract 


Flow of viscous fluids are not usually discussed in detail in general and basic 
courses of physics. This is due in part to the fact that the Navier-Stokes 
equation has analytical solution only for a few restricted cases, while more 
sophisticated problems can only be solved by numerical methods. In this text, 
we present a computer simulation of wind tunnel, f.e., we present a set of 
programs to solve the Navier-Stokes equation for an arbitrary object inserted 
in a wind tunnel. The tunnel enables us to visualize the formation of vortices 
behind object, the so-called von Karman vortices, and calculate the drag force 
on the object. We believe that this numerical wind tunnel can support the 
teacher and allow a more elaborate discussion of viscous flow. The potential 
of the tunnel is exemplihed by the study of the drag on a simplihed model of 
wing whose angle of attack can be controlled. A link to download the programs 
that make up the tunnel appears at the end. 


Keywords: fluid dynamics, viscous flows, wind tunnels, drag force. 


Introducao 

/ 

Pode uma bolinha em queda frear? Um experimento recente (1) mostra-nos 
que, surpreendentemente, a resposta e sim. Os autores sugerem que isto se deve a 
formagao gradativa de uma esteira de vortices no ar atras da bolinha. Estes vortices, 
que surgem em virtude da viscosidade do fluido, sao arrastados pela bolinha em 
queda e, por isso, sao responsaveis pela forga de arrasto que se opoe ao movimento. 
No inicio do processo de formagao da esteira, a forga de arrasto e proporcional a ve- 
locidade. Uma vez que a esteira tenha sido completada, a forga se torna proporcional 
ao quadrado da velocidade. Na transigao entre um regime e outro, o arrasto torna-se 
maior que o peso e a frenagem e observada na bolinha. Este cenario e conhrmado 
por meio de simulagoes numericas (2, 3). 

Embora interessantes e proficuos, veja por exemplo (4, 5), problemas fisicos 
que envolvem dinamica de fluidos viscosos nao sao, em geral, discutidos em cursos de 
fisica basica e geral. Quando isso ocorre, o assunto e tratado de maneira superheial 
(6, 7). Por que? Para responder a esta pergunta precisamos primeiro dehnir o que 
sao fluidos viscosos e discutir como sao descritos. 

Em mecanica dos fluidos, costuma-se dehnir como huido qualquer substancia 
que evidentemente nao seja um solido. Quando submetido a uma tensao de cisalha- 
mento (tipo de tensao gerada por forgas que atuam em sentidos opostos) um huido 
tende a se deformar. A viscosidade, uma propriedade do huido que determina quao 
facilmente o huido se deforma, e uma fungao do huido em questao e da temperatura 
e corresponde ao atrito interno resultante das interagoes entre as moleculas (7, 8). 
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Fluidos sao descritos por meio de uma equagao derivada independentemente 
por G. G. Stokes (1816 - 1903) e L. M. H. Navier (1785-1836), a equagao de Navier- 
Stokes (9) k 

^ = ^v=a-vx(n><.-), (1) 

em que F e o campo de velocidades, G = V x F e a vorticidade e Re e o numero de 
Reynolds, uma grandeza adimensional resultante da divisao de forgas de inercia por 
forgas de viscosidade em um elemento de fluido. Explicitamente, 


Re 


pvD 

R 


( 2 ) 


em que v e a velocidade do vento (longe do obstaculo, para o caso de um objeto 
inserido no interior de um tiinel de vento), D e uma dimensao caracteristica do ob¬ 
jeto, /i e a viscosidade e p e a massa especifica do fluido. O numero ou coeflciente de 
Reynolds foi introduzido por G. G. Stokes em 1851. Gontudo, foi O. Reynolds (1842 
- 1912) que popularizou o seu uso em 1883 e demonstrou, pela primeira vez, que 
o numero poderia ser utilizado como um criterio para distingao entre escoamento 
laminar, no qual as trajetorias das particulas do fluido tendem a ser paralelas, e tur- 
bulento, no qual as trajetorias das particulas sao curvas, irregulares e entrecruzadas 
(8, 10, 11). Por exemplo, para o escoamento sobre um cilindro, a transigao entre o 
regime laminar e o turbulento ocorre quando 140 < Re < 300 (12, 13). 

A equagao de Navier-Stokes admite solugoes analiticas para alguns poucos 
casos. Na verdade, a prova matematica da existencia de uma solugao global para a 
equagao de Navier-Stokes ainda nao existe e e um dos Millennium Prize Problems 
(7). Gom isso, resta-nos recorrer as simulagoes numericas. Embora metodos numeri- 
cos para solugao da equagao de Navier-Stokes sejam numerosos na literatura (veja, 
por exemplo, (14, 15, 16, 17)), simulagoes voltadas para o ensino sao bem menos 
abundantes. Neste texto apresentamos um tiinel de vento simulado no qual se pode 
escolher o formato do objeto a ser inserido no tiinel assim como a velocidade do vento. 
Imaginamos que esta ferramenta seja capaz de subsidiar o professor e permitir que 
uma discussao mais profunda de escoamentos viscosos seja realizada. 

Este texto esta organizado da seguinte maneira. Inicialmente, os procedimen- 
tos para solugao da equagao de Navier-Stokes para um objeto de formato arbitrario 
no interior de um tiinel de vento sao apresentados. Em seguida, descrevemos um 
metodo, recentemente introduzido (3), para determinagao da forga sobre objetos 
imersos em tiineis de vento. Posteriormente, exempliflcamos potencialidades didati- 
cas colocando no interior do tiinel de vento uma asa cujo angulo de ataque pode ser 
controlado. Finalmente, em uma segao concludente, apresentamos alguns comenta- 
rios que julgamos ser relevantes ao leitor e eventual utilizador do tiinel simulado. 

^ Pode-se encontrar facilmente na literatura dedugoes fenomenologicas da equagao de Navier- 
Stokes, embora nao haja uma dedugao a partir de primeiros principios aplicados a estrutura 
atomica microscopica do fluido. Veja, por exemplo, (7, 9). 
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1 Solucao da equacao de Navier-Stokes para urn objeto de 
formato arbitrario imerso em um tunel de vento 

O sistema fisico em questao e um tunel de vento, no qual a velocidade do 
vento e controlavel e onde pode ser inserido um objeto de formato arbitrario, por 
exemplo, uma asa on um cilindro. 

Para descrever a dinamica deste sistema fisico, resolvemos a equagao 1 em 
uma rede bidimensional de 400 x 200 pontos (correspondente a uma regiao que deno- 
minamos regiao de interesse), suficiente para nossos propositos (a resolugao podera 
ser ampliada no caso de necessidade de maior precisao). Para resolver a equagao 
1, usamos um metodo (2) de diferengas finitas com relaxagoes sucessivas sobre as 
vorticidades As configuragoes iniciais dos campos de velocidades e vorticidades 
sao as resultantes da solugao da equagao 1 no limite em que Re —> 0. Neste limite, a 
equagao de Navier-Stokes reduz-se a equagao de Laplace. Utilizar esta configuragao 
inicial e necessario pois do contrario, um comportamento transiente espurio pode 
ser obtido. Veja, por exemplo, (22). 

As condigoes de contorno sao: (a) na superficie do obstaculo e dentro dele, 
sao mantidas nulas as velocidades e vorticidades; (b) fora da regiao de interesse 
(definida acima), o vento tern apenas uma componente nao nula, ao longo da diregao 
X, de modo que nesta regiao, v = Vx (por conveniencia, fez-se V = 1 em unidades 
adimensionais) e as vorticidades sao nulas. 

Em principio, objetos de qualquer formato podem ser inseridos no tunel res- 
guardada a condigao de que suas dimensoes nao ultrapassem 20% do tamanho total 
do tiinel, por exemplo, ou fragao menor ainda caso seja necessario maior precisao. 
Ademais, o objeto precisa estar localizado no centro do tiinel. Atendidas estas con¬ 
digoes, fomos bem sucedidos em observar a formagao da esteira de vortices para Re 
< 1000, sendo possivel simular com niimeros de Reynolds proximos de 1000 sem 
que alteragoes significativas no tamanho da rede, nas condigoes de contorno e na 
viscosidade numerica sejam necessarias ^ . A figura 1 mostra um vortice surgindo 
atras de um objeto imerso no tiinel, neste caso, uma asa. No link: <https://www. 

^ Aplicamos o metodo das relaxagoes sucessivas aqui para resolver as equagoes de Stokes {Re = 0) 
e de Navier-Stokes {Re 7 ^ 0). No primeiro caso, o metodo consiste simplesmente na substituigao 
da vorticidade em um ponto da rede pela media das vorticidades dos seus quatro vizinhos 
imediatos. Para o caso da equagao de Navier-Stokes, outros termos aparecem e estes tambem 
precisam ser expresses em fungao dos valores das velocidades e vorticidades nos pontos em 
seu entorno, tantos pontos quantos forem necessaries para que todas as diferengas finitas que 
substituirao as derivadas tenham precisao pelo menos ate a segunda ordem no passo discrete 
espacial e no passo de tempo. Para uma explanagao do uso de relaxagoes sucessivas na solugao 
de equagoes diferenciais, veja (18, 19, 20). Curiosamente, na referenda (21, Problem 3.29), o 
metodo das relaxagoes e aplicado na solugao de um problema de eletrodinamica, o que traz 
a tona o conhecido fato de que a hidrodinamica e a eletrodinamica sao formalmente muito 
semelhantes. 

^ A introdugao da viscosidade numerica e um procedimento muito usado tradicionalmente para 
suavizar as variagoes espaciais bruscas no campo de velocidades. Trata-se simplesmente de 
acrescentar as vorticidades uma pequena fragao de seu proprio Laplaciano. A viscosidade nu¬ 
merica foi utilizada por nos na solugao da equagao de Navier-Stokes para Re > 0. Sua utilizagao 
produz uma suavizagao numerica, mas pode introduzir desvios quantitativos no resultado. Em 
principio, a fragao do Laplaciano somada ao campo de vorticidades deve ser pequena o sufi- 
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Figura 1 - Vortice aparecendo atras de um objeto imerso no tiinel. 


dropbox.eom/s/ls4jbk589dfl4qp/tunel asa 1000 22.wmv?dl=0> pode-se baixar 
um video onde e mostrada a evolugao do campo de velocidades quando, partindo 
da configuragao de Stokes, o vento e ligado em Re = 1000 e escoa sobre uma asa 
com o angulo de ataque acionado em 22°. A evolugao temporal observada atesta a 
adequagao de nossos resultados a outros, experimentais e computacionais, presentes 
na literatura (9, 12, 15, 23, 24). 

Contudo, como e possivel determinar a forga de arrasto sobre o objeto no 
tiinel? A proxima segao e dedicada a descrigao de um metodo geral capaz de faze-lo 
(3). 

2 Um metodo alternativo para determinacao da forca sobre 
objetos imersos em tuneis de vento 

A maneira tradicional de calcular a forga de arrasto sobre um obstaculo ri- 
gido envolve o conhecimento de dois elementos basicos: o campo de pressoes e a 
viscosidade nas proximidades da interface objeto/fluido, be., na regiao correspon- 
dente a camada limite (25). E possivel, alternativamente, obter a forga de arrasto 
unicamente a partir do conhecimento do campo de velocidades do fluido em torno 
de um obstaculo, por meio do gradiente deste campo ao longo da superficie da inter¬ 
face objeto/fluido. Neste caso, deve-se realizar uma integragao sobre esta superficie 
(8, 26, 27). Quando o campo de velocidades e determinado em pontos de uma rede 
discreta, a precisao para o calculo do gradiente fica comprometida, a menos que se 
adote uma rede muito fina perto da superficie de integragao, o que requer um grande 
esforgo computacional. Para contornar este problema tecnico, propomos substituir 

ciente para que tais desvios estejam controlados. No caso de Re = 1000, a fragao adotada foi 

1 / 100 . 
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a integral de superficie por uma integral de volume sobre todo o volume ocupado 
pelo objeto, como descrevemos a seguir. 

No instante t, o campo de velocidades em cada ponto r da grade e ni(r). 
Esta configuragao no instante t e obtida a partir de configuragoes em instantes an- 
teriores t — dt, t — 2dt, etc. Nos pontos da grade no interior do obstaculo e na sua 
superficie, Ut = 0 em qualquer instante t. Supontiamos que o objeto seja removido 
no instante t e o volume que o objeto ocupava seja preenchido com fluido estatico. 
Partindo da configuragao do campo de velocidades ja conhecido Vt{r), obtem-se a 
configuragao futura no tempo t + Dt . Como o objeto rigido foi substituido por um 
fluido, algumas velocidades nao nulas aparecem no interior do volume anteriormente 
ocupado pelo objeto. Em outras palavras, a remogao do objeto permite que o campo 
de velocidades penetre um pouco no interior do sen volume. Agora, este campo de 
velocidades interno pode ser integrado no volume. O resultado e multiplicado pela 
densidade do fluido, obtendo-se assim o momento que seria transferido do fluido 
para o objeto. Dividindo-se esse impulse por Dt, finalmente obtem-se a forga de 
arrasto Pode-se observar na figura 2 a forga de arrasto sobre um longo cilindro 
estatico perpendicular ao vento no interior do tiinel, em fungao do tempo; esta forga 
corresponde aproximadamente a situagao fisica que descrevemos na segao introdu- 
toria deste texto. O comportamento da forga de arrasto para o cilindro estatico se 
harmoniza com a interpretagao sugerida na segao introdutoria para o problema da 
bolinha em queda: ha um tempo transiente durante o qual a esteira de von Karman 
ainda nao se formou e a forga de arrasto nao e proporcional a o que ocorre tao 
logo a esteira esteja formada. Na transigao entre estes regimes, a forga de arrasto 
torna-se maior que o peso da bolinha e a mesma experimenta uma frenagem. 

Naturalmente, o metodo acima descrito pode ser aplicado nao apenas a um 
cilindro, mas a um objeto de formato qualquer. 


3 Uma asa em um tunel de vento 

Para exemplihcar as potencialidades do tunel, apresentamos, nesta segao, um 
estudo do comportamento da forga de arrasto sobre um modelo simplihcado de asa 
cujo angulo de ataque, i.e., o angulo que o eixo da asa faz com o vento, pode ser 
controlado. 

Nossa simulagao corresponde a situagao fisica de um aeromodelo, cuja asa 
tern largura caracteristica aproximada de 60cm voando a uma velocidade aproxi- 
mada de 24m/s no ar, cuja viscosidade e de 1.4207.10“®m^/s. A asa que desenha- 
mos e mostrada na hgura 3. Para a velocidade do vento hxa em Re = 1000, o 
angulo de ataque pode variar entre ±22° sem que as paredes do tunel interhram 
fortemente na simulagao. Neste caso, apos o regime transiente, vortices girando em 
sentidos contraries aparecem alternadamente atras da asa, assim como ocorre com 
o cilindro. 

^ O tempo Dt necessario para que o vento entre na regiao anteriormente ocupada pelo objeto 
varia de acordo com o mimero de Reynolds e com as condigoes iniciais. Considera-se que uma 
estimativa do tempo e aceitavel, quando multiplicado por dois, o momento transferido torna-se 
duas vezes maior. 
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Figura 2 - Forga de arrasto sobre um cilindro estatico no tiinel de vento. Inicial- 
mente o tiinel esta desligado. Em seguida, em t = 0, o tiinel e ligado 
com niimero de Reynolds Re = 1000. A componente da forga na diregao 
perpendicnlar ao flnxo do vento e desprezivel. A forga de arrasto e repre- 
sentada em nnidades arbitrarias, ignorando fatores como a densidade do 
flnido, etc. Dnrante nma nnidade de tempo, o vento viaja o eqnivalente 
ao diametro do cilindro. A componente da forga qne e paralela ao flnxo 
de vento anmenta, atinge nm valor maximo, diminni e finalmente se es- 
tabiliza, o qne corresponde a nma sitnagao dinamica em qne aparecem 
continnamente vortices sncessivos girando em sentidos alternados atras 
do cilindro. A esteira de vortices de von Karman e entao formada. O 
destaque mostra o mesmo nnma escala longa no tempo. 



Figura 3 - Modelo simplificado de asa introduzida no tiinel. O angulo de ataque 
pode ser controlado entre ±22°. 


A componente longitudinal do arrasto neste caso apresenta um comporta- 
mento qualitativo analogo ao observado para o cilindro (veja figura 2): a forga cresce, 
atinge um valor maximo, reduz-se um pouco e estabiliza, passando a oscilar em torno 
de um valor fixo positivo. Este valor fixo anmenta com o angulo de ataque. O corn- 
portamento da componente transversal do arrasto, igualmente, depende do angulo 
de ataque. Para um angulo de ataque igual a 0°, observa-se qne esta componente, 
depois de um estagio transiente, flutua continnamente em torno de zero, o qne e 
esperado para um obstaculo axialmente simetrico (8). A medida qne o angulo de 
ataque cresce, os valores desta componente aumentam. Este resultado corresponde 
satisfatoriamente ao qne se sabe sobre a relagao entre a sustentagao, forga corres- 
pondente a componente transversal do arrasto, e o angulo de ataque (28, 29). Veja 
as figuras 4 e 5. 

Em nossa simulagao, a asa pode ser ligeiramente modificada para incluir um 
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Figura 4 - Componentes do arrasto para angulo de ataque igual a 0°. Em (a), 
apresenta-se a componente longitudinal e, em (b), a componente transver¬ 
sal. A componente longitudinal cresce, atinge um valor maximo, reduz-se 
um pouco e se estabiliza, passando a oscilar levemente (imperceptivel na 
figura) em torno de um valor fixo positivo, enquanto a componente trans¬ 
versal, apos o transiente, passa a fiutuar em torno de um valor fixo, no 
caso do angulo de ataque igual a 0°, este valor e nulo, como esperado 
para obstaculos simetricos. 



Figura 5 - Componentes do arrasto para angulo de ataque igual a 15°. Em (a), 
apresenta-se a componente longitudinal e, em (b), a componente trans¬ 
versal. Assim como antes, a componente longitudinal cresce, atinge um 
valor maximo, reduz-se um pouco e estabiliza, passando a oscilar em 
torno de um valor fixo positivo enquanto a componente transversal, apos 
o transiente (neste caso, bem diferente do caso para angulo de ataque 
0°, pois a simetria axial foi quebrada), passa a fiutuar em torno de um 
valor fixo, no caso do angulo de ataque igual a 15°, este valor e positivo, 
embora pequeno se comparado com a componente longitudinal. 
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Figura 6 - Modelo simplificado de asa com flap. 



Tempo 

Figura 7 - Comparagao das componentes do arrasto para a asa com e sem flap para 
um angulo de ataque igual 15°. Em (a) e (b), se observam, respectiva- 
mente, as componentes longitudinal e transversal do arrasto para asa 
sem flap (repetido da flgura 5) enquanto em (c) e (d), se observam, res- 
pectivamente, as componentes longitudinal e transversal do arrasto para 
asa com flap. Como esperado, para um mesmo angulo de ataque, a asa 
com flap apresenta valores mais altos para suas componentes embora o 
comportamento qualitative seja aproximadamente o mesmo com on sem 
flap. 


tipo de flap. Flaps consistem em abas on superficies articuladas existentes na parte 
posterior das asas que alteram temporariamente a geometria das mesmas. Quando 
abaixados e/ou estendidos, os flaps aumentam o arrasto e a sustentagao, o que e 
especialmente importante no procedimento de pouso. O flap que instalamos na asa 
e mostrado na flgura 6 e e equivalente a um flap do tipo ventral. As alteragoes nas 
componentes do arrasto em fungao da presenga do flap para um angulo de ataque 
igual a 15° sao mostradas na flgura 7. Da analise do graflco, percebe-se que a presenga 
do flap altera a forga de arrasto, aumentando o valor de ambas as componentes em 
relagao ao sen valor sem o flap, o que se harmoniza com o que se sabe a respeito da 
fungao deste tipo de flap nas asas (8). 


Consideracoes finals 

/ 

O estudo de escoamentos viscosos, embora riquissimo, e praticamente ausente 
em cursos de fisica geral e basica. Isso se deve, em parte, ao fato deste tipo de 
escoamento ser descrito pela equagao de Navier-Stokes que, por sua vez, so admite 
solugao analitica para alguns poucos casos. Diante deste quadro, apresentamos neste 
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texto um tiinel de vento simulado no qual se pode escolher o formato do objeto a 
ser inserido no tiinel assim como a velocidade do vento. O objetivo do tiinel e 
instrumentalizar o professor e viabilizar que uma discussao mais aprofundada do 
tema seja realizada. 

Os programas que compoem a simulagao, assim como um arquivo com ins- 
trugoes para utilizagao do tiinel podem ser baixados no enderego <https://www. 
dropbox. com / s/t3y8og20zfliypk8 / programas. rar? dl=0 > . 

Alguns comentarios finals se fazem necessaries. O tiinel, concebido para ana- 
lisar escoamentos bidimensionais, pode ser facilmente modificado para dar conta de 
escoamentos tridimensionais. Porem, nossa experiencia com a utilizagao do tiinel 
indica que isso exigiria enorme esforgo computacional, o que ate entao nao pode ser 
contornado nem mesmo com a utilizagao de clusters de alto desempenho. Isso se 
deve ao fato de ser dificil realizar um processo de computagao em paralelo para o 
tiinel uma vez que uma configuragao e gerada a partir da anterior. 

No tiinel, “as paredes” que limitam a regiao de interesse, inevitavelmente, 
interferem nos resultados. Para minimizar sua infiuencia, seria necessario adotar 
paredes muito mais afastadas do objeto. Na verdade, para Re pequenos, teriamos 
adicionalmente que adotar um tiinel tridimensional, porque em duas dimensoes, a 
infiuencia das paredes se faz sentir ate distancias proporcionais a l/i?e (30, 31). Esta 
interferencia, contudo, nao modifica os principals aspectos qualitativos da solugao, 
aspectos estes que sao o foco deste trabalho. 

Num formato piloto, o tiinel foi utilizado em duas disciplinas, fisica termica 
para o curso tecnico de automagao industrial e fisica geral II para licenciatura em 
fisica do Instituto Federal do Rio de Janeiro. Na aplicagao no curso de automagao, 
por tratar-se de um curso de nivel medio, o topico escoamentos viscosos foi apre- 
sentado como complementar. Neste caso, o problema da bolinha que freia ao cair 
foi apresentado e discutido em detalhes com o auxilio do tiinel. Simulagoes com 
asas tambem foram utilizadas para discutir as razoes fisicas pelas quais o aviao 
consegue manter-se no ar. No curso de licenciatura, o tiinel foi utilizado por um 
pequeno niimero de alunos para elaboragao do projeto de final do curso. Neste caso, 
inicialmente, os alunos se familiarizaram com tiinel numerico. Posteriormente, inse- 
riram diferentes objetos no tiinel, como uma bala de revolver, um triangulo e um 
cilindro girando com velocidade angular controlavel. Os dados obtidos por eles para 
este ultimo caso foram contrapostos com os presentes em (3), o que viabilizou uma 
interessante discussao sobre o efeito Magnus e suas manifestagoes. 
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